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We investigate the response of a circumstellar accretion disc to the fly-by of a 
perturbing mass on a parabolic orbit. The energy and angular momentum transferred 
during the encounter are calculated using a reduced three-body method. In almost 
all close encounters the energy and angular momentum transfer is dominated by disc 
material becoming unbound from the system, with the contributions from close disc 
particle - star encounters being significant. For more distant encounters with some 
prograde element to the motion the disc material loses energy and angular momentum 
to the perturber's orbit through a resonance feature. The magnitude of the energy 
transfer calculated in our simulations is greater than that of the binding energy of 
material exterior to periastron by a factor of two in the prograde case, and up to a 
factor of five in the case of the retrograde encounter. The destructive nature of the 
encounters indicates that a non-linear treatment is essential in all but the most distant 
encounters. 

Key words: accretion, accretion discs - celestial mechanics, stellar dynamics - binaries: 
close - stars: formation. 
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2 Much of the recent work on binary star formation has 
stressed the importance of interactions between stars and 
their associated gaseous protostellar discs. A proposed mech- 
anism for the formation of binary and multiple systems 
is that of 'prompt initial fragmentation' (Pringle 1989) in 
which stars form in bound groups, each comprising a hand- 
ful of protostellar objects at separations of the order a Jeans 
length. Each fragment then undergoes the processes of sin- 
gle star formation to yield a central star surrounded by a 
protostellar disc. At the same time the stars interact un- 
der their mutual gravitation. The weak clustering seen af- 
ter a few million years in star-forming regions is consistent 
with such a formation mechanism. Optical observations of 
Taurus-Auriga (Gomez et al. 1993) reveal six such clusters 
containing approximately 15 stellar systems in each. Speckle 
imaging surveys of the same region (Ghez et al. 1993) sug- 
gest that a large fraction ( ^ 60 per cent) of T Tauri stars 
have companions at projected separations comparable with 
typical disc radii, implying that star-disc interactions may 
have been important in their formation. In this paper we 
present a new method for calculating the energy and an- 
gular momentum transfer in such star-disc encounters and 
discuss the implications of our results for binary star forma- 
tion. 



It should be noted that collision rates determined 
from number densities and velocity dispersions observed in 
present-day clusters are rather low, and we must therefore 
appeal to an earlier epoch in the cluster's formation his- 
tory when number densities were higher, velocity dispersions 
were smaller, and collision rates were substantial. Based on 
the prompt formation scenario outlined above, one would ex- 
pect stellar number densities to be around a few x 10^ pc~^ 
and velocity dispersions in the sub-1 km regime. Also, 
the mass of a protostellar disc is much larger during the 
earliest stages of evolution since much of the enshrouding 
material has yet to fall in, and thus it is at these earliest 
stages that discs have the greatest effect on the orbital evo- 
lution of the system. The presence of a massive disc offers 
an enhanced collisional cross-section to nearby members and 
hence a mechanism for the transfer of orbital energy and an- 
gular momentum between perturber orbit and disc material. 
Penetrating encounters also offer an explanation for the en- 
hanced weak-line T- Tauri binary distribution at small sepa- 
rations (Ghez et al. 1993) by removing much of the material 
at large radii. 

Recent estimates of the rate of binary star formation 
through star-disc interactions have modelled such encoun- 
ters in a prescriptive fashion (Clarke & Pringle 1991a, 1991b; 
Murray, Clarke & Pringle 1991; Larson 1990). For exam- 
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pie, Clarke & Pringle (1991a) took the energy transferred 
to be equal to the binding energy of the disc material ex- 
terior to periastron, and showed, in large- systems, that 
even the densest protostellar regions could only yield a bi- 
nary fraction of up to 20 per cent by this method. In subse- 
quent papers (Clarke & Pringle 1991b; McDonald & Clarke 
1995) the authors followed the evolution of a small number 
of protostars and their discs, experimenting with a number 
of prescriptions similar to that described above. They found 
that in small- clusters the presence of discs assists purely 
dynamical capture and can give rise to a substantial binary 
fraction. The binary statistics were, however, somewhat sen- 
sitive to the prescription employed, thus underlining the ne- 
cessity of obtaining a reliable prescription for this process. 

In the case of distant encounters, analytical results for 
energy and angular momentum transfer can be used, these 
being based on linear perturbation theory. All such theories 
rely on the perturbing mass to excite waves and resonances 
in the disc material and calculate the energy and angular 
momentum transfer induced as a result. Such analysis (Os- 
triker 1994; Papaloizou & Terquem, in preparation) is simi- 
lar in both method and results to that on tidal interactions 
of galaxies. Ostriker calculates energy and angular momen- 
tum transfer between a protostellar disc and a perturbing 
star on a parabolic orbit over a range of disc radius to pe- 
riastron separation ratios and disc-orbit orientations. She 
finds that for close encounters disc angular momentum is 
reduced slightly, and that capture will only occur for po- 
lar and retrograde encounters. Papaloizou & Terquem use a 
similar analysis in computing the angular momentum flux 
between a disc and a perturbing body and find that the disc 
loses angular momentum to the perturber and thus accretion 
on to the central star is increased, and that for non-coplanar 
encounters the evolution of the system is not necessarily to- 
wards coplanarity. 

For close and penetrating encounters linear theory is no 
longer appropriate and different calculational techniques are 
required. Limited numerical simulations of penetrating star- 
disc encounters (Clarke & Pringle 1993, hereafter CP93) 
have shown that outside periastron most of the disc ma- 
terial becomes unbound and thus a suitable estimate of the 
energy transferred might be the binding energy of the disc 
material outside periastron. Some attempt was made to ac- 
count for viscosity using a 'sticky-particle' formulation and 
it was found that the principal effect of the viscosity was that 
it enabled re-circularization of the remaining disc material 
after the fly-by. 

Similar studies were conducted by Heller (1993) using 
Smoothed Particle Hydrodynamics (SPH), where the effect 
on the stellar orbit was also calculated. To date, this rep- 
resents the only such calculation in the non-linear regime. 
This study did not, however, explore the question of which 
area of the disc is chiefly responsible for transfer of energy 
and angular momentum to the orbit, and therefore raises 
the question of whether the critical regions of the disc were 
adequately resolved. A further uncertainty concerns the in- 
exactitude of the hydrodynamic forces in SPH on a particle 
by particle basis, a problem that is negligible in many ap- 
plications where the bulk properties of the flow are of inter- 
est, but which may be important in the present problem, if 
such effects shift particles in and out of critical regions of 
the disc. On the other hand, a Lagrangian method, such as 



SPH, is essential in such calculations, as grid-based methods 
soon founder, during penetrating disc encounters, due to the 
enormous dynamic range in densities that develop in tidally 
disturbed regions (Korycansky & Papaloizou 1995). 

Faced with these difficulties, and given also the current 
interest in obtaining reliable star-disc interaction prescrip- 
tions, we have opted for an approach that begins by dis- 
secting the star-disc interaction process in the absence of 
hydrodynamic forces. The main motivation of this work is 
to discover the critical regions of the disc that couple to the 
stellar orbit and to classify the particle orbits that dominate 
the energy and angular momentum transfer. Such a study 
can be used as groundwork for future hydrodynamical cal- 
culations, by indicating the regions of the disc that need 
to be most closely resolved. Furthermore, there are reasons 
for believing that the results obtained in this ballistic case 
may not be drastically different from those that would apply 
when hydrodynamical effects are included: the simulations 
of CP93 show that the majority of particles in the outer 
disc are stripped into diverging orbits by the encounter, so 
that ballistic modelling of the flow is not likely to be greatly 
erroneous in this case. We therefore, in this paper, examine 
the ballistic case, since it is clearly analy sable on a parti- 
cle by particle basis and because it permits the economical 
method (described in Section 2). We bear in mind, however, 
that our present conclusions may be somewhat modified by 
gas dynamical calculations. 

Our treatment of star-disc interactions in this paper re- 
calls early investigations of interacting galaxies (Toomre & 
Toomre 1972; Wright 1972; Fneev, Kozlov & Sunyaev 1973; 
Yabushita 1977) although, due to the computing time re- 
quired, these studies were mainly qualitative investigations 
emphasizing the morphology of distorted galaxies. In the 
present paper we use a new method for calculating the en- 
ergy and angular momentum transferred in a star-disc en- 
counter, based on the reduced three-body problem. Since we 
are ignoring all inter-particle interactions in this model, as 
well as the effect of the disc material on the binary orbit dur- 
ing the numerical experiment, particles may be integrated 
one at a time. At the end of the experiment we calculate the 
energy and angular momentum changes of all the particles 
and equate this to the corresponding change in the binary 
orbit. Section 2 lays out the basis for this approach, whilst 
Section 3 outlines the numerical method and Section 4 de- 
scribes our results. In Section 5 we apply our results to some 
typical systems, and we draw our conclusions in Section 6. 



2 DERIVATION OF REDUCED THREE-BODY 
RESULTS 

Consider a system of three masses M2 , Mi and m (a typical 
disc particle) located at R2, Ri and R and moving with 
velocities V2, Vi and V with respect to the centre of mass 
of the system. We may then write the total angular momen- 
tum, J, of the system as 

J = M2(R2 X V2) Mi(Ri X Vi) m(R X V) (1) 
and the total energy, E, as 



iM2V2 + -MiVl 

2^2 ^ 



GM2M1 , 1 ^^2 

— — — - H — mV 

R2 - Ri 2 
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GMim 
R,i — R, 



GM2m 
R-2 — R 



(2) 



We need to know the changes in the energy and angu- 
lar momentum of the binary orbit as it evolves under the 
influence of the surrounding disc material. Since we ignore 
the effect of the disc material on the binary orbit, we must 
evaluate the changes in the energy and angular momentum 
of the disc material itself. These may then be equated, with 
the appropriate change of sign, to the change in energy and 
angular momentum of the binary orbit, since total energy 
and angular momentum must be conserved. If this is cal- 
culated directly, however, we run into problems. Consider 
the case when a disc particle, m, is bound to one of the 
stars, say Mi, while the other star, M2, recedes. The angu- 
lar momentum of the particle in the centre of mass frame 
is given by m(R x V). As IR-I^- 00, V tends to the orbital 
velocity of m about Mi which is finite, and, for a general 
non-circular orbit, periodically variable in both magnitude 
and direction. Thus, as M2 recedes, the disc particle's angu- 
lar momentum oscillates with a frequency equal to that of 
the m — Ml system and grows in amplitude. Since the disc 
particle's angular momentum is clearly time-dependent, we 
cannot guarantee a correct value for this (or other similar 
quantities) at the end of the simulation. A method which 
removes this inconvenience will now be described. 

Initially, all disc particles are substantially closer to Mi 
than to M2 (assuming a suitable starting time has been cho- 
sen) whilst, after the fly-by, each disc particle may be in one 
of three situations. 

(i) Mass m remains substantially closer to Mi than 
to M2. 

(ii) Mass m is substantially closer to M2 than Mi . 
(iii) Mass m is substantially further from Mi and M2 
than they are from each other. 

In each case we can split the three-body system into two 
'two-body' subsystems - the close pair of masses becomes 
one system and its centre of mass and the third mass become 
the other. By writing the energy and angular momentum in 
a way appropriate to each of the three possible situations we 
are able to obtain correct values for the energy and angular 
momentum of the system. 

2.1 Case (i): mass m remains substantially closer 
to Ml than to M2 

We first define some convenient coordinates. Let rc and Vc 
be the position and velocity of the centre of mass of the 
m — Ml system defined by 



(Ml -I- m)rc = MiRi -I- mR, 
(Ml -I- m)vc = MiVi -I- mV, 



(3a) 
(3b) 



and ri, r and Vi, v be the positions and velocities of Mi and 
m in the centre of mass of the m — Mi system; thus 



Ri = rc-l-ri, 
R = rc -I- r, 

Vi = Vc -I- Vi, 
V = Vc+V. 

Combining (3) and (4) gives the obvious 



(4a) 
(4b) 
(4c) 
(4d) 



Miri -I- mr = 0, 
MiVi +mv = 0. 

Equation (1) now becomes 



(5a) 
(5b) 



J« = M2(R2 X V2)-FMi[(rc-Fri) X (vc-Fvi)] 

-I- m [(rc -I- r) X (vc -|- v)] . (6) 

Rearranging and noting (5) yield 

J«=J«+J« (7) 

where 

J*'' = M2(R2 xV2) + {Mi+ m)(rc x Vc), (8a) 

jW = Mi(ri X vi) -I- m(r X v). (8b) 

Jq' is the angular momentum of M2 and of the centre 
of mass of the Mi — m system treated as a point mass in the 
centre of mass frame. 

Jj'' is the angular momentum of Mi and m about their 
mutual centre of mass. Note that these two terms are pro- 
portional to Mim^ and Mj^m respectively and thus in the 
limit that m <^ Mi only the latter need be considered. 

Equation (2) now becomes 



fiW = ^M2Yl + -Ml [vc - 



2 



2 



Ml 



GM1M2 



rc - R2 + ri 



1 r Ml 1^ GMim 
+ 0"^ vi r 

2 L m J I ri — r I 

GM2m 



I rc - R2 + r I 

Rearrangement yields 



(9) 



rnfi) 1,^,,2 1/,^ X 2 1 m^V^ 1 Mi^Vi 

E^ ' = -M2V2 -I- -(Ml -I- m)Vc -I — 1 

2 ^^2^ '^" 2Mi 2 m 



GM1M2 



GMim 

rc - R2 + ri I I ri - r I 

GM2m 



R2 +r 



(10) 



We now use the condition that Mi and m are much 
closer to each other than either is from M2 , or, more strictly, 
that |r|, |ri| ^ |rc — R2I. Thus we can expand the fifth and 
seventh terms to give: 



where 



=iM2V2' + i(Mi+m)vc' 



2 GM2(Mi-Fm) 



R2 



^(i)_lm^v^ IMi^v? GMim 
^ ~ 2 Ml 2~ I ri - r I' 



(11) 



(12a) 



(12b) 



Eg^ is the total energy of M2 and of the Mi — m system 
treated as a point mass. is the 'internal' energy of the 
Ml — m system. The first two terms in (12b) are proportional 
to Ml and mMf respectively and thus in the limit that 
m ^ Ml only the latter need be considered. 
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2.2 Case (ii): mass m is substantially closer to M2 
than to Ml 

The second case is identical to the first but for an inter- 
change of subscripts, thus: 

= , (13) 

where 



J<"' = Mi(Ri X Vi) + (M2 + m)(rc x Vc 

J*"' = M2(r2 X V2) + m(r x v), 

and 



here 

Bf' = iMiV? + i(M2 + m)v, 



2 GMi{M2+m) 



rc - Ri 



(ii) _ 1 m^v^ lM|vi GM2m 



2 M2 2 m I r2 - r 



(14a) 
(14b) 

(15) 

(16a) 
(16b) 



Note that the definitions of Vc,rc,van(ir change ac- 
cordingly and that V2 and r2 have the obvious meaning. 
Also, the appropriate terms may be omitted provided that 
m < M2. 



2.3 Case (iii): mass m is substantially further from 
Ml and M2 than they are from each other 

Case (iii) is a little different. 

Define rc and Vc to be the position and velocity of the 
centre of mass of the Mi — M2 system defined by 



(Ml + M2)rc = MiRi M2R2, 
(Ml + M2)vc = Ml Vi + M2V2, 



(17a) 
(17b) 



and ri,r2 and Vi,V2 to be the positions and velocities of 
Ml and M2 in the centre of mass of the Mi — M2 system; 
thus 



Ri = rc-l-ri, 

R2 = rc -I- r2 , 

Vi = Vc -I- Vi, 

V2 = Vc -I- V2 . 



(18a) 
(18b) 
(18c) 
(18d) 



Combining (17) and (18) gives the obvious 
Min -FM2r2 = 0, (19a) 
Mivi -I- M2V2 = 0. (19b) 

Equation (1) now becomes 



= M2[(rc + r2) x(vc+V2)] 

+ Ml [(rc ri) X (vc vi)] -Fm(R X V) 

Rearranging and noting (19) yield 

j(m) ^ j(m) 

where 

jf =M2(r2 X V2)-FMi(riX,vi), 

jf = (Ml + M2)(rc X Vc) + m(R x .V). 



(20) 

(21) 

(22a) 
(22b) 



Jq"' is the total angular momentum of the Mi — M2 sys- 
tem in its centre of mass frame, j''"' is the angular momen- 
tum of m and the Mi — M2 system treated as a point mass in 
the centre of mass frame. The two terms in (22b) are propor- 
tional to Mm^ and M^m respectively (where M = Mi -I-M2) 
and in the limit that m <^ M only the latter need be con- 
sidered. 

Equation (2) now becomes 



1 ,^ / ^2 1 ,^ / x2 GM1M2 
7-M2 Vc -FV2r + 77M1 Vc Vlf - -. '-^ 

2 ' 2 ' ri - r2 



2 



GMim 



GM2m 



rc -I- ri - R I I rc -I- r2 - R 

Rearrangement yields 

= ^(Mi + M2)wl + iM2vi + ^Miv? 



(23) 



GMim 



GM1M2 , 1 ^,2 

1 r H — V — — - 

I ri - r2 I 2 I rc ri - R 

GM2m 



(24) 



I rc r2 - R 

We now use the condition that Mi and M2 are much 
closer to each other than either is to m, or, more strictly, 
that |ri|, |r2| ^ |rc — R|. Thus we can expand the sixth and 
seventh terms to give: 



where 



-MiVf -I- -M2V 



2 



GM1M2 
I ri - r2 I 



(25) 



(26a) 



^r' = >+M2)vg+imv-- ^(;^-tgf • (26b) 

i?o"'' is the total energy of the binary orbit in its centre 
of mass frame. is the total energy of m and of the Mi — 

M2 system treated as a point mass. The first two terms in 
(26b) are proportional to Mm^ and mM^ respectively and 
in the limit that m <^ M only the latter need be considered. 



In the limit m ^ Mi we may take Bq'' ~ E^ 



^0 

and Jq ' ~ Jq"' ~ Jq i so therefore the changes in energy 
and angular momentum of the system are simply equal to 
the difference between the final and initial values of Ei and 
Ji respectively, independent of which approach is required to 
calculate the initial and final values. These changes may then 
be directly equated to the changes in the energy and angular 
momentum of the binary orbit, from which the final binary 
orbital parameters may be obtained from the standard 

GM^i 



2a 



J = n^jGMa(l - e2). 



(27a) 



(27b) 



where M is the total mass, M = Mi -|- M2, fj. is the reduced 
mass of the binary system, jj, = — jj-^, a is the semi-major 
axis and e is the eccentricity. 



3 NUMERICAL METHOD 

The disc particles initially describe Keplerian orbits about 
the primary. Mi . Since the particles are non-interacting, the 
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results are independent of the number of disc particles used, 
as long as the mass attributed to any one particle is much 
less than either of the stellar masses. In this limit, the en- 
ergy and angular momentum changes of a disc particle are 
proportional to its mass, thus allowing any disc density pro- 
file to be fitted later, provided that the aforementioned in- 
equality holds. We must also be sure that there are enough 
particles throughout the disc to follow any localized regions 
of interest. Disc particles are massless and thus the response 
of the binary to the disc is not calculated during the ex- 
periment. No form of dissipative process is employed and 
therefore the disc particles are integrated ballistically in the 
time-dependent binary potential. The binary consists of two 
point masses, the primary. Mi, and the perturber, M2, on 
parabolic orbits about their centre of mass. We begin the 
calculation at a time of t = to (t = at periastron) when 
the binary separation is at least ten times the initial disc 
radius, fdisc. The initial positions are found by solving Ke- 
pler's equations (Taff 1985). We find that the disc evolution 
is insensitive to the initial binary separation as long as the 
initial separation is substantially larger than the initial disc 
radius. 

The motion of both the disc particles and the point 
masses is integrated using the Bulirsch-Stoer method. We 
use an unsoftened gravitational potential for each star. The 
accuracy of this method was checked by integrating the sys- 
tem over the time of a typical fly-by, but using a zero-mass 
perturber. We found energy conserved to better than one 
part in 10^ and angular momentum to one part in 10^^. 

At the end of each timestep we determine into which 
of the three categories described above each of the particles 
falls, and Ei and Ji are calculated appropriately (equations 
8b, 12b, 14b, 16b, 22b, and 26b). Only when Ei and Ji have 
stabilized and we can be sure that our values for the energy 
and angular momentum transfer are sufficiently accurate do 
we stop integration of that particle. Particles are integrated 
in a serial fashion. 

Particles that lie within some chosen radius of Mi or M2 
and are gravitationally bound to that mass are considered 
accreted to that mass and are no longer integrated. Since 
our chosen accretion radius of 0.5 per cent times the initial 
disc radius, i.e. faccrete = 5 xl0~^rdi3c, is always much less 
than our interstellar separation we can be sure that the val- 
ues calculated for Ei and Ji at the point where we decide 
the particle has been accreted will be accurate. Thus, even 
though particles may be accreted early in the simulation, we 
can confidently include their contribution to the transfer. 



4 RESULTS 

We present the results of our tests of the numerical method 
and of several models in which we vary the disc radius with 
respect to fperi and the relative inclination of the disc and 
orbital planes. In all models M2/M1 = 1 and the binary 
orbit is parabolic. Our code time units are 

r = i(r^,„/GM)^/^ 

= O.O8(rpe„/au)'/'(M/M0)-'/' yr, (28) 

where M = Mi + M2. 

Figure 1 shows the axes we use. The primary lies at 



Figure 1. We employ axes in which the primary lies at the origin 
and the perturber orbits in the x — y plane. The binary angular 
momentum vector, Jorbit 5 is therefore in the positive z direction. 
The disc's angular momentum vector, Jdisc 5 lies in the y — z plane 
and makes an angle, (?, with that of the orbit in the sense shown. 

the origin and the perturber orbits in the x — y plane. The 
binary angular momentum vector is therefore in the positive 
z direction. The angular momentum vector of the disc, Jdisc, 
always lies in the y — z plane and makes an angle, 9, with 
that of the binary, Jorbit, in the sense indicated. ^ = 0° in 
the prograde encounter and 180° in the retrograde scenario. 

4.1 The numerical method 

The accuracy of the Bulirsch-Stoer method has already been 
established and we do not go into the details here. In or- 
der to ensure that we have enough particles in our disc, we 
choose a representative ring undergoing a strong interaction 
with the perturber and monitor the change in its energy and 
angular momentum as we vary A^, the number of particles 
that it contains. We use a ring of radius l.Qfperi and sev- 
eral values for N, ranging from 10 to 1000. At low values 
of N, accuracy is very poor, especially in the angular mo- 
mentum, whilst for N > few x 10^ both energy and angular 
momentum are within a few per cent of the values obtained 
with 1000 particles in the ring, an indication that we are 
populating the disc densely enough to be sure that we are 
following all of the interaction. Unfortunately, the compu- 
tational time required to evaluate 1000 particles in this ring 
implies prohibitively long computational times for an entire 
disc, and thus we choose N = 100 as a compromise between 
computational speed and numerical accuracy. Having said 
that, however, when N ~ 100, the tolerance in the energy 
is less than 1 per cent and less than 10 per cent in angular 
momentum. We have also examined rings at other radii and 
found that, for similar accuracy, 50 particles are sufficient 
at r/fperi = 1.0 and a mere 13 in the innermost region at 
r/fperi = 0.25. 

4.2 Coplanar, prograde encounter 
4-2.1 Particle distributions 

In the same vein as CP93 we first discuss the results for 
a coplanar, prograde encounter. The stars are on parabolic 
orbits and M2/M1 = 1. Our post-encounter particle dis- 
tributions are similar to those of CP93 and the fractions 
of particles that become unbound or are captured (Fig. 2) 
are also in good agreement with CP93. The section of Fig. 
2 for r/fperi < 0.8 should be compared to fig. 3 of CP93. 
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Figure 2. The fate of particles as a function of initial radius 
for a coplanar, parabolic, prograde orbit with the mass ratio 
M2/M1 = 1. At each radius within the disc, the fraction of parti- 
cles at that radius that remain bound to the primary (asterisks), 
that are captured by the perturber (filled squares) and that be- 
come unbound from the system (open squares) are shown. 



Figure 3. The fraction of disc particles in each of the three final 
states, bound to the primary (asterisks), captured (filled squares) 
and unbound (open squares), plotted against time, t; t = at 
periastron. The encounter is coplanar, parabolic and prograde 
with M2/M1 = 1. 



Outside periastron the number of particles that are cap- 
tured falls to only a few per cent compared with 40 per 
cent at r/fperi = 0.6. The amount of material remaining 
bound to the primary reaches a minimum of 20 per cent at 
r/rperi = 1.3 which rises again to 40 per cent at the disc's 
edge. The remainder of the material becomes unbound, the 
exact amount varying between 60 and 75 per cent, peaking 
at r/fperi = 1.5. The fact that this fraction is so high adds 
weight to the validity of the assumption made in previous 
papers (McDonald & Clarke 1995; Clarke & Pringle 1991b) 
that the reaction is dominated by the unbinding of all the 
disc material exterior to periastron. 

Fig. 3 shows the time evolution of these fractions. The 
most important thing to notice is that the bulk of the inter- 
action takes place within 1 code unit of periastron (t = 0), 
when the perturber is well within rdisc of the primary, and 
by a time of t = 5.0 the final state of the particles is evi- 
dent. We find that 58 per cent of the disc particles become 
unbound, 30 per cent remain bound to the primary and the 
remaining 12 per cent are captured by the perturber. These 
percentages reflect the number of test particles in our disc 
(using a uniform surface number density) and not disc mass 
which depends on the disc's surface density profile through 
the radial dependence of the mass ascribed to each particle. 

4-2.2 Energy and angular momentum transfer 

Since we find that the majority of particles have settled into 
their final bound or unbound states very soon after perias- 
tron passage, further integration is only necessary to ensure 



that we have accurate values for the energy and angular 
momentum according to Section 2. We find that, by a time 
t = 300, over 90 per cent of particles have been integrated 
for a sufficient length of time to ensure that their energies 
and angular momenta have settled to within the tolerance 
level specified (one part in 10^). The remaining 10 per cent 
take a little longer to complete, although not prohibitively 
so. Fig. 4 shows the particles in their initial positions with 
respect to the primary at {0,0}. The magnitudes of the en- 
ergy and angular momentum changes are represented by the 
area of the symbol. 

We find, as we might expect, some structure in the pat- 
tern of energy transfer over the disc's area. In the central 
regions (r/fperi ~ 0.3) there is very little change, whilst at 
intermediate radii (0.3 ^ r/fperi ~ 1.5) the energy transfer is 
fairly uniformly distributed. At large radii (r/fperi ~ 1.5) we 
find there is a significant excess confined to a small range in 
azimuthal angle, this range being different at different radii. 
It is clear that most of the interaction is not taking place 
in one small region of the disc and thus there is no need to 
overpopulate parts of the disc selectively. 

The peaks in the energy transfer as shown in Fig. 4(a), 
form a spiral pattern and one might expect these peaks to 
become aligned at some point due to the effects of differ- 
ential rotation associated with the Keplerian nature of the 
disc. However, the picture is somewhat more complicated 
than that. Selecting the dominant particles only (Figure 5a), 
and following their motions as they approach periastron, we 
found that they form four distinct groups (Fig. 5b). Group 1 
is composed of approximately 20 per cent of these particles, 
which undergo close encounters with the perturber as their 
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Figure 4. The locations of disc particles in the disc at the start 
of the simulation. The area of each symbol is proportional to the 
absolute magnitude of (a — top) the energy and (b — bottom) 
the angular momentum changes that the particle suffers during 
the encounter. The encounter is coplanar, parabolic and prograde 
with M2/M1 = 1. 

orbital motion about the primary brings them close to the 
perturber (Fig. 6a). Group 2 particles, again comprising ap- 
proximately 20 per cent, undergo close encounters with the 
perturber as its motion brings it close to them (Fig. 6b). 
The third is by far the dominant group, comprising 50 per 
cent of particles. As the perturber approaches the disc, these 
particles are disrupted from their original orbits and go into 
orbit about the perturber in the opposite sense to that of the 



Figure 5. (a — top) The peaks in energy transfer, selected from 
Fig. 4(a). The system is shown at time t = —40. (b — bottom) 
The same system at a time t = —2, showing how the particles 
have evolved into four distinct groups. Groups are labelled left to 
right — group 1 (open triangles), group 2 (filled triangles), group 3 
(open squares) and group 4 (filled squares). The perturber's orbit 
(indicated by the solid line) is prograde (anti-clockwise here), 
parabolic and coplanar. The central star and the perturber are 
represented by the bold, open squares. The area of the symbol is 
proportional to the magnitude of that particle's energy change. 
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Figure 6. The time evolution of tlie pealis in energy transfer sliown in Fig. 5. (a — top left) At time t = —1.6. (b — top right) At time 
t = -1.1. (c - bottom left) At time t = -0.5. (d - bottom right) At time t = -0.3. 



disc about the primary. They then undergo close encoun- 
ters with the primary and become unbound (Fig. 6c). The 
remaining particles form the fourth group and, like those 
in the third group, are sufficiently perturbed that they be- 
gin to orbit the primary in the opposite sense and suffer 
close encounters with the perturber as it approaches perias- 
tron (Fig. 6d). Note that in Figs 6(c) and (d) the particles 
from groups 1 and 2 can be seen scattered throughout the 
system and their motion away from the centre of mass of 
the system is evident. These peaks carry approximately one 
third of the total energy transferred, suggesting that close 
encounters of the kind demonstrated here are significant and 
warrant closer examination. It is clear that this is a highly 
non-linear process and thus must be treated numerically. 

There is no obvious correlation between peaks in en- 
ergy transfer and peaks in angular momentum transfer 
throughout the disc, except in the almost undisturbed, cen- 
tral (r/fperi ~ 0.3) region. It is therefore imperative that 



we follow all of the disc throughout the interaction, rather 
than just concentrate on the energy or angular momentum 
'hotspots'. Binning particles according to the magnitude of 
the change in their energy and angular momenta (Fig. 7), we 
see that the majority of the particles have relatively small 
changes and that these particles produce the dominant ef- 
fect (Fig. 8). This is reassuring in as much has already 
been said, there is no one small region of the disc, and thus 
only a few particles, that dominates the interaction. 

We can see clearly that interior to r/fperi ~ 0.3 there is 
almost no perturbation to the disc material. Fig. 9(a) shows 
the energy change at each radius. Between r/fperi ~ 0.3 and 
Si 0.6 the energy change is negative; the material between 
these radii is becoming more tightly bound and, correspond- 
ingly, the binary is becoming more weakly bound. Thus, if 
our disc were not to extend farther than 0.6rperi (this value 
depending on the mass ratio M2 /Mi , and the initial binary 
orbital parameters), we would expect the binary system to 
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Figure 7. The particles are binned according to tlie magnitude of 
tlieir energy and angular momentum changes during a coplanar, 
prograde, parabolic encounter with M2/M1 = 1. The numbers of 
particles in each bin are plotted. The solid line indicates binning 
by energy, the dashed line by angular momentum. 



gain energy. Note that co-rotation between the perturber 
and disc material occurs at r/fperi = 0.63 at periastron and 
we suspect that the negative-going part of Fig. 9(a) may 
be a co-rotation resonance feature. It has been pointed out, 
however, that an m = 2 inner Lindblad resonance (occur- 
ring at r/fperi = 0.4) is also possible (E.G. Ostriker, private 
communication). We are unable to distinguish between the 
two with our current results. It is interesting to note that the 
inner part of this feature is dominated by material that re- 
mains bound to the primary - at this radius approximately 
95 per cent of the material remains bound - whilst the outer 
part of the resonance is dominated by material that is cap- 
tured by the perturbing mass. At this radius (r/fperi ~ 0.5) 
approximately one third of the material is captured (Fig. 
2). Thus even when rperi/^disc ~ 2 the interaction is already 
non-linear and can no longer be treated in a linear fashion. 

Beyond r/fperi ~ 0.6, the curve rises sharply to a peak 
just inside periastron. In this region, as one might expect 
given the destructive nature of the interaction, the reaction 
is now dominated by the material that becomes unbound 
from the system. There is some residual effect from material 
that remains bound to the system, acting in the opposite di- 
rection, i.e. it loses energy to the binary orbit. In the region 
1 ^ r/fperi ~ 2 the total energy transferred is approximately 
constant, coinciding with a peak in the fraction of mate- 
rial becoming unbound from the system (Fig. 2). Beyond 
r/fpeii ~ 2, the total energy drops as approximately r~^, 
being almost wholly due to the unbound material; material 
that ends up being bound to either star contributes very lit- 
tle to the total energy transfer. The fluctuations present in 
the curve in this region have been shown to be \/N noise. At 



Figure 8. The particles are binned according to the magnitude of 
their energy and angular momentum changes during a coplanar, 
prograde, parabolic encounter with M2/M1 = 1. The total contri- 
butions to the energy and angular momentum transfer from the 
particles in each bin are plotted. The solid line indicates binning 
by energy, the dashed line by angular momentum. 

these radii, the dependence of the energy is reflected in 
the fraction of material becoming unbound (Fig. 2). These 
results agree with those of Yabushita (1977) who found that 
material at radii r/fperi = 0.375 and r/fperi = 0.5 acceler- 
ated the relative motion of the perturber whilst material at 
r/rperi = 0.625 tended to decelerate the relative motion of 
M2 and Ml . 

Fig. 9(b) shows the equivalent results for angular mo- 
mentum. Again within r/fperi ~ 0.3 the transfer is negligi- 
ble. The same resonance feature is present between r/fperi ~ 
0.3 and r/fperi ~ 0.6 where the disc material is losing an- 
gular momentum to the binary orbit - the structure of this 
feature is as for the energy transfer. Outside r/fperi ~ 0.6 
the angular momentum transfer rises sharply, again domi- 
nated by the unbound material. The total transfer is approx- 
imately constant out to r/fperi ~ 1.5 where it falls equally 
sharply to zero at r/fperi ~ 2. Here, the transfer is equally 
dominated by the unbound material gaining angular mo- 
mentum and the bound material losing angular momentum. 
Beyond r/fperi ~ 2, the total angular momentum transfer 
becomes increasingly negative with increasing r, only it is 
now dominated by the material that remains bound. This is 
coincidental with an increasing fraction of disc material re- 
maining bound at these radii (Fig. 2) and reinforces our ear- 
lier findings showing no correlation between peaks in energy 
and angular momentum transfer. The angular momentum 
transfer at large radii is dominated by material that suffers a 
reduction in its angular momentum whilst remaining bound. 
Its contribution to the energy transfer at these radii is small, 
however, due to its small initial energy. There is a 'peak' in 
the total angular momentum transfer at r/fperi ~ 2.8 which 



10 S.M. Hall, C.J. Clarke and .I.E. Pringle 



Figure 10. The fate of particles as a function of initial radius 
for a coplanar, parabolic, retrograde orbit with the mass ratio 
M2/M1 = 1. At each radius within the disc, the fractions of par- 
ticles at that radius that remain bound to the primary (asterisks), 
that are captured by the perturber (filled squares) and that be- 
come unbound from the system (open squares) are shown. 



The fractions falling into the three final configurations are as 
follows: unbound 26 per cent, bound to the primary 52 per 
cent and captured by the secondary 22 per cent. Similarly 
the energy transfers for each of the three cases are tend- 
ing towards zero. In the case of the angular momentum, 
things are not quite so clear. The contributions from pri- 
mary and secondary bound material remain at values sim- 
ilar to those at r/fperi = 4 whilst the transfer due to the 
unbound material changes sign although the magnitude re- 
mains approximately the same. The consequence of this is 
that the total angular momentum transfer maintains an ap- 
proximately linear increase in magnitude to a value of —1.35 
at r/fperi = 10. 



Figure 9. The energy change (a — top) and angular momentum 
change (b — bottom) per unit mass plotted against radius for the 
coplanar, parabolic, prograde encounter with M2/M1 = 1. The 
contributions from primary-bound material (asterisks), captured 
material (filled squares) and unbound material (open squares) are 
plotted along with the total (open triangles). 

is also present in the transfer from the material bound to 
the primary. There is no evidence that this is a resonance 
feature. Captured material is also important at large radii, 
its effect to some extent cancelling that of the unbound ma- 
terial. 

It would appear that by r/fperi ~ 3 the transfer curves 
are tending to the limiting values at large radius. Calcula- 
tions with a ring of material at r/fperi = 10 confirm this. 



4.3 Coplanar, retrograde encounter 

In this encounter we now have the disc rotating in the op- 
posite sense to the perturber about the primary star. In this 
case we not expect to find a resonance feature, for obvious 
reasons. 

4-3.1 Particle distributions 

Again, our post-encounter particles distributions are in good 
agreement with those of CP93. Overall (using a constant 
surface density profile), we find that 62 per cent of the 
disc particles become unbound from the system during the 
encounter, 28 per cent remain bound to the primary and 
the remaining 10 per cent are captured by the perturbing 
body. Fig. 10 shows how the fractions of particles in each 
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Figure 11. The locations of disc particles in the disc at the start 
of the simulation. The area of each symbol is proportional to the 
absolute magnitude of (a — top) the energy and (b — bottom) the 
angular momentum changes that the particle suffers during the 
encounter. The encounter is coplanar, parabolic and retrograde 
with M2/M1 = 1. 

of the three final states vary with the initial radius. Within 
r/fperi ~ 0.7 all disc particles remain bound to the pri- 
mary. Outside this, the number that become unbound rises 
sharply, peaking at 75 per cent at r/fperi ~ 1-6. Correspond- 
ingly, the fraction that remain bound drops sharply, reach- 
ing a minimum of 25 per cent at the same radius. This fig- 
ure should be compared with fig. 8 of CP93 for the region 
r/fperi < 2. Beyond r/fperi ~ 2 some material is captured. 



Figure 12. (a — top) The peaks in energy transfer, selected from 
Fig. 11(a). The system is shown at time t = —40. (b — bottom) 
The same system at a time t = —2, showing how the particles 
have evolved into three distinct groups. Groups are labelled 1 to 
3 in a clockwise fashion starting from the left — group 1 (open 
triangles), group 2 (filled triangles) and group 3 (open squares). 
The perturber's orbit (indicated by the solid line) is parabolic and 
coplanar. Disc material orbits in a retrograde fashion (clockwise 
here). The central star and the perturber are represented by the 
bold, open squares. The area of the symbol is proportional to the 
magnitude of that particle's energy change. 
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Figure 13. The time evolution of tlie pealis in energy transfer sliown in Fig. 12. (a — top left) At time t = —1.2. (b — top right) At time 
t = 0.5. (c - bottom) At time t = 0.8. 



the fraction increasing linearly with radius, reaching 25 per 
cent at r/fperi = 4. The fraction remaining bound to the 
primary drops with increasing radius in this region, reaching 
20 per cent at r/fperi = 4. Correspondingly, the fraction 
becoming unbound remains approximately constant exterior 
to r/fperi ~ 2 at approximately 65 per cent. 

4-3.2 Energy and angular momentum transfer 

Fig. 11 shows how the energy and angular momentum trans- 
fers are distributed throughout the material. Again, there is 
a noticeable pattern in the energy transfer: the inner, undis- 
turbed region is larger than before, extending to r/fperi ~ 
0.8. In the region 0.8 ^ r/fperi ~ 2 there is a large amount 
of energy transfer whilst at larger radii the peaks in energy 
transfer are localized within small ranges of azimuthal an- 
gle, different at different radii. These peaks form a spiral 
structure as was present in the prograde case, but of oppo- 



site sense here. This again suggests an alignment of these 
peaks via the differential rotation of the disc (clockwise in 
the retrograde case). 

Selecting these peaks (Fig. 12a) we find that they form 
not one but three distinct clumps by a time t = —2 (Fig. 
12b). The third group is the most dominant of the three, con- 
sisting of approximately 40 per cent of the most dominant 
particles. Their motion around the central star brings them 
close to the perturber as it traverses the disc where they un- 
dergo close encounters with the perturber, being scattered in 
the process. It should be noted that, due to the proximity of 
their encounters, none is scattered through angles less then 
7r/2 (Fig. 13a). The second group comprises around half of 
the particles which deviate from their orbits around the pri- 
mary due to the presence of the perturber and undergo close 
encounters with the primary (Fig. 13b). The first group is 
much smaller, comprising only 10 per cent of the particles 
which are scattered by the perturber as it passes through 
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the disc on its way away from the central star - again there 
are no scatterings through an angle less than n/2. 

As for the prograde case, we find no correlation between 
peaks in energy and angular momentum transfer, except for 
the relatively undisturbed inner region. The angular mo- 
mentum transfer is spread approximately uniformly over the 
region exterior to r/fperi = 2. 

Fig. 14(a) shows the energy transfer as a function of 
radius for the coplanar, retrograde encounter. The undis- 
turbed region near the primary now extends to r/fperi ~ 0.8. 
We expect this since the relative velocity of the particles 
with respect to the perturber in this region is high, so the 
time available for the interaction is short as well as the 
separation being large. We see no resonance feature, as ex- 
pected. This undisturbed region confirms the calculations 
of Yabushita (1977) who found energy changes caused by 
material interior to r/fperi = 0.625 to be insignificant. In 
the region 0.6 ^ r/fperi ~ 1 there is a small blip in the en- 
ergy transfer performed by the primary-bound material. In 
this region, however, the energy transfer rises very steeply 
to a peak at r/fperi ~ 1.2; the energy transfer at this maxi- 
mum being over twice that at the maximum in the coplanar 
case. The energy transfer is approximately constant over the 
range 1.2 ^ r/fperi ~ 1.6. At large radii the energy transfer 
falls off as approximately . Throughout, the total energy 
transfer is overwhelmingly dominated by the unbound ma- 
terial. Even though an increasing fraction is being captured 
at large radii (Fig. 10), its effect is negligible throughout the 
range examined. 

Fig. 14(b) shows the equivalent result for the angular 
momentum transfer. Note that the initial angular momen- 
tum of the binary is positive, that of the disc is negative, and 
so the angular momentum of the disc material is increasing 
in the direction of that of the orbit. Again, the inner region 
(»"/rperi ~ 0.8) is undisturbed. At r/fperi ~ 0.8, the angular 
momentum transfer rises steeply and is dominated by un- 
bound material. The rise becomes less steep at periastron 
and continues approximately linearly with radius towards 
large radii. Outside periastron the transfer is still dominated 
by the unbound material although not as completely as at 
periastron. Fxterior to r/fperi ~ 1.6 there is a significant 
contribution from the material remaining bound to the cen- 
tral star, remaining approximately constant in this region. 
At large radii (r/fperi ~ 3), there is an increasing contribu- 
tion from captured material, corresponding to the increased 
fraction of material being captured at these radii. 

4.4 Angle of encounter 

While, for a binary system forming out of a single rotating 
cloud or from disc fragmentation, one might expect there to 
be a correlation between the angular momenta of the disc 
and of the binary orbit, we may not get such a correlation 
in binaries formed during the collapse of a cluster of stars. 
Surveys of solar-type binary systems reveal similar correla- 
tions (Hale 1994, for example) which show that for binary 
systems there is a correlation between the equatorial plane 
and the orbital plane only for a ^ 30 au and that there is no 
preferred orientation in hierarchical systems. It is therefore 
important that non-coplanar encounters are considered. 

We now consider three different relative orientations of 
the disc and orbit, 9 = 45°, 90°, and 135°, where prograde 



Figure 14. The energy change (a — top) and angular momentum 
change (b — bottom) per unit mass plotted against radius for the 
coplanar, parabolic, retrograde encounter with M2/M1 = 1. The 
contributions from primary-bound material (asterisks), captured 
material (filled squares) and unbound material (open squares) are 
plotted along with the total (open triangles). 

encounter is at ^ = 0° and retrograde is at ^ = 180°. In all 
cases the binary's major axis lies in the disc plane (Fig. 1). 

We recognize the importance of encounters where pe- 
riastron does not lie in the disc plane, but due to space 
restrictions we reserve these for a follow-up paper. However, 
following some preliminary calculations, we do not expect a 
substantial difference between orientations with periastron 
outside the disc plane and those presented here. 
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Figure 15. Energy (a — top left) and the three components [x — (h — top right), y — (c — bottom left), ^ — (d — bottom right)) of the 
angular momentum transfer per unit mass in the three non-coplanar encounters: 45° open triangles, 90° open squares and 135° open 
circles. The angular momentum of the original binary orbit is in the positive z direction. 



4.4.1 Particle distributions 



4.4.2 Energy transfer 



Using a constant surface density profile in all three cases we 
find the fractions of disc particles that end up in the three 
possible situations are as follows. In the 9 = 45° encounter, 
we find that, overall, approximately 33 per cent of the disc 
particles become unbound from the system, 32 per cent re- 
main bound to the primary and 35 per cent are captured by 
the perturber. For 9 = 90°, 57 per cent of the disc particles 
become unbound from the system during the encounter, 36 
per cent remain bound to the primary and 7 per cent are 
captured. These fractions are somewhat similar to those for 
the 9 = 135° case at 53, 46 and 1 per cent respectively. 



In all three encounter orientations we can see that, as be- 
fore, the inner regions are only slightly disturbed (Fig. 15). 
For 9 = 45°, a resonance feature is present although at a 
slightly different radius from the prograde case. This may 
be expected since we must now consider the projection of 
the motion of the disc material on to the plane of the bi- 
nary orbit. For 9 = 45° co-rotation occurs at r/fperi = 0.5 
and the m = 2 inner Lindblad resonance at r/fperi = 0.3. 
There is a small negative-going region at r/fperi = 0.5 in the 
9 = 90° encounter, probably caused by some warping of the 
inner disc plane during fly-by that enables a small resonance 
effect to occur. 

Fxterior to r/fperi ~ 0.6, the energy transferred rises 
sharply to a peak exterior to periastron, in the range 
1 ^ r/fperi ~ 2, before falling off at larger radii. The 
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Figure 16. The total (a — top) energy change and (b — bottom) 
angular momentum change of the disc material within a radius r 
per unit mass within a radius r plotted against r for the coplanar, 
prograde encounter. The binary orbit is parabolic and M2/M1 = 
1. Each of the three density profiles is shown: constant surface 
density (open triangles), cx (open squares), cx r~^l^ (open 

circles). The filled symbols indicate the binding energy of the 
material between that radius and r/rperi = 1 for the respective 
profiles. 



Figure 17. The total (a — top) energy change and (b — bottom) 
angular momentum change of the disc material within a radius r 
per unit mass within a radius r plotted against r for the coplanar, 
retrograde encounter. The binary orbit is parabolic and M2 /Mi = 
1. Each of the three density profiles is shown: constant surface 
density (open triangles), cx (open squares), cx r~^l^ (open 

circles). The filled symbols indicate the binding energy of the 
material between that radius and r/rperi = 1 for the respective 
profiles. 
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magnitude of the peak in energy transfer is similar for all 
three encounter angles and comparable to that in the pro- 
grade case. In all cases the energy transfer is dominated by 
the material that becomes unbound from the system. 

4.4-3 Angular momentum transfer 

In the 9 = 45° encounter, the angular momentum in the x 
direction (along which the binary's major axis lies) hardly 
changes (Fig. 15b). The transfer is dominated in the outer 
regions (r/fperi ~ 2.4) by material whose y component is 
being reduced (Figure 15c). (The initial disc angular mo- 
mentum is in the y + z direction.) Interior to r/fperi ~ 1.5 
the transfer is dominated by material whose z component is 
increasing (Fig. 15d). In between, both effects are approxi- 
mately equal. 

In the 9 = 90° orientation the disc's initial angular 
momentum is in the y direction. From Figs 15(b)-(d) it 
is clear that there is little change in any component in the 
innermost regions of the disc whilst, at intermediate radii 
(0.8 ^ r/fperi ~ 2), the transfer is dominated by the z com- 
ponent which is becoming more positive at the expense of 
the binary's angular momentum (originally in the positive 
z direction). At large radii (r/fperi ~ 2), we find that the 
y component peaks strongly at r/fperi ~ 3.3 (Fig. 15c) and 
dominates the interaction. The z component undergoes a 
small change in this region whilst the change in the x com- 
ponent is of intermediate value. 

We find that all three components are of approximately 
equal magnitude at large radii in the 9 = 135° encounter 
whilst, in the inner regions, changes in the x and y com- 
ponents dominate. The initial disc angular momentum is in 
the y — z direction. 

In all cases material gains angular momentum in the 
z direction and thus the z component of the initial binary 
orbital angular momentum is reduced. As was found in the 
coplanar encounters, the interaction is dominated by mate- 
rial becoming unbound from the system. 



5 APPLICATION OF RESULTS 

The quantitative results we have presented so far have been 
'per unit mass', thus allowing any reasonable radial density 
profile to be assumed later. In order to fit an arbitrary disc 
density profile, S(r), it is simply a matter of multiplying the 
radial energy and angular momentum transfer curves by 

E(r)27rrAr (29) 

where Ar is the radial extent of each annulus of material. 
We now apply our results to some 'typical' discs. It should 
be noted that we still scale our masses by Mi and distances 
by fperi and any application to real systems will be left as 
an exercise for the reader. One caveat to this is that in the 
regime of large disc masses, i.e. Mdisc/Mi ^ 1, this linear 
relationship no longer applies, not only because the dynam- 
ics of the disc would be affected by its self-gravity, but also 
because the effect of the disc material on the binary orbit 
would be significant, thus invalidating the reduced three- 
body approach (Section 2). 

We present the results of three density profiles in such 
a way as to simplify their application to specific systems. A 



constant surface density was chosen as reference profile and 
contrasted with the results for two power-law profiles. 

Figs 16(a) and (b) show the energy and angular momen- 
tum changes for the prograde encounter for each of the four 
profiles. Plotted vertically is the energy change of the mate- 
rial interior to radius r divided by the mass in the same re- 
gion. Thus in order to apply a specific disc mass, Mdisc / Mi , 
and radius, r/fperi, one simply reads the energy or angular 
momentum change per unit mass off the appropriate profile 
curve at the required value of r/fperi and multiplies by the 
ratio Mdisc /Mi . Also plotted is the binding energy of the 
material exterior to periastron as given by 

/''disc QJ^ 
E 27rrdr. (30) 
2 r 
peri 

Figs 17(a) and (b) show the equivalent results for the retro- 
grade encounter. 

It should be noted that, for the steeper density profiles, 
the magnitude of the energy change is smaller than for other 
profiles since a steeper profile results in more of the mass 
being present in the inner regions where the transfer effects 
are smaller. 

In both prograde and retrograde cases and for all den- 
sity profiles, the energy transferred is always greater than 
the binding energy of the disc. The energy transfer is a fac- 
tor of approximately two larger than the energy transfer in 
the prograde case, rising to a factor of five in the retrograde 
case. Thus any assumption that equates the energy loss of a 
binary orbit to the binding energy of disc material exterior to 
periastron is underestimating the strength of the encounter. 
This effectively means that the disc masses in previous pre- 
scriptive work (Clarke & Pringle 1991a, 1991b; Murray et 
al. 1991, for example) can be reduced by a factor of between 
two and five. We stress, however, that the results obtained 
here for low-mass discs may not necessarily be applicable to 
the massive discs (Mdisc ~ Mi) considered in these works. 

Interior to periastron, the energy and angular momen- 
tum changes are far from small and when inside r/fperi ~ 
0.6, where the resonance feature dominates, the energy and 
angular momentum of the disc material are actually reduced. 
It is therefore imperative that the effects of material interior 
to periastron are accounted for, particularly if periastron lies 
near the edge of the disc. 

Tables 1 and 2 show the results of applying Figs 16 and 
17 to some 'typical' disc systems in the prograde and retro- 
grade cases respectively. We choose two typical disc masses, 
Mdisc /Ml = 10"^ and Mdisc /Mi = 10~\ and disc radii, 
»"disc/»"peri = 0.5 and rdisc/»"peri = 2. In all cases we have 
calculated the parameters of the final binary orbit as per 
equations (27). The initial orbital energy and angular mo- 
mentum take the values Binitiai = (implying a parabolic 
orbit) and Jinitiai = 0.5z (indicating periastron separation) 
respectively. 

In all cases, we find that the periastron separation of 
the final binary orbit is very close to that of the original. 
The maximum deviation shown in our examples is of 14 
per cent for a constant-E disc with mass Mdisc /Mi =0.1 
and rdisc/»"peri = 2 in the retrograde case. This, combined 
with the abrupt nature of the interaction, suggests that an 
impulse approximation applied at periastron would not be 
drastically erroneous. 
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Table 1. The application of our results to the three disc density profiles, flat, cx r ^ and cx r ^/'^ . We use two disc masses and two disc 
radii. The encounter is coplanar and prograde with M2/M1 = 1. Our initial energy and angular momentum take the values .^'initial = 



and Jinitial = 0-5 respectively, a' is the semi-major axis, e' is the eccentricity and 1 
binary orbit. 
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Table 2. The application of our results to the three disc density profiles, flat, cx r — 1 and cx r— . We use two disc masses and two disc 
radii. The encounter is coplanar and retrograde with M2/M1 = 1. Our initial energy and angular momentum take the values il/initial = 
and Jinitial = 0-5 respectively, a' is the semi-major axis, e' is the eccentricity and r'^^^-^ is the periastron separation of the post-encounter 
binary orbit. 



S(r) 


"•disc 
Tperi 


Mdi 


-bit 


Mdi 


■bit 




g' 




Ml 

e' 


/ 

r 

pen 
Tperi 


g' 

Tper 


Mdi 
M 


, = c ^ 10-1 

1 

e' 


/ 

r 

pen 
Tperi 


const. 


0.5 


-2.86 X 


: 10-3 


-8.17 X 


: 10" 


-4 


1.75 X 


10" 


0.999943 


0.999988 


1.75 X 


103 


0.99943 


0.999960 




2 


-1.31 : 


X 10" 


-6.35 X 


: 10" 


-1 


3.81 X 


lOl 


0.974 


0.988 


3.81 X 


10" 


0.775 


0.859 


cx r — 1 


0.5 


-3.84 X 


: 10-3 


-4.48 X 


: 10- 


-4 


1.30 X 


104 


0.999923 


1.000020 


1.30 X 


103 


0.99923 


1.00020 




2 


-9.10 X 


: 10-1 


-4.13 X 


: 10" 


-1 


5.49 X 


lOl 


0.982 


0.9925 


5.49 X 


10" 


0.833 


0.919 


—7/4 

cx r ' 


0.5 


-8.04 X 


: 10-3 


-1.26 X 


: 10- 


-4 


6.22 X 


103 


0.99984 


1.000076 


6.22 X 


102 


0.9984 


1.00075 




2 


-3.13 X 


: 10-1 


-1.33 X 


: 10" 


-1 


1.60 X 


102 


0.9938 


0.9978 


1.60 X 


lOl 


0.939 


0.977 



The sign of the change in periastron of the orbit may 
be predicted from ABorbit and AJorbit using the standard 
equations of orbital mechanics to obtain 
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for our parabolic orbit. For an increase in periastron, i.e 
Afperi > 0, then we have: 

f > i if ABorbit > 
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AJorbit 
A-Borbit 



< i if ABorbit > 
> i if ABorbit < 



Fig. 18 shows the ratio AJorbit /ABorbit of the material 
at radius r plotted against r for both the prograde and ret- 
rograde encounters. In both cases the inner material affects 
the orbit in both senses. In the prograde case, material be- 
tween r/fperi = 0.6 and r/fperi = 1.8 acts so as to decrease 
fperi, just as One might expect since the orbit is becom- 
ing more tightly bound. Outside r/fperi = 1.8, however, the 



material acts so as to increase fperi. This may be a little 
surprising but is due to the large increase in orbital angular 
momentum. In the retrograde case all material exterior to 



r/fpe 



0.8 acts so as to decrease fperi as one might expect. 



We have not plotted AJorbit /ABorbit at the smallest radii as 
we begin to suffer numerical error. It should be noted that 
in these regions (i.e. r/fperi < 0.3) linear theory becomes 
appropriate. 

Manipulation of the standard orbital equations also 
yields: 



ABorbit AJorbit 



2eAe 



Jot 



(1- 



(34) 



For initial orbits that are close to being parabolic, i.e. 
e ~ 1 or i? ~ 0, then the AJorbit /^orbit term is small com- 
pared with the ABorbit /Borbit term and can be ignored, so 
that 



Ae oc AEo 



(35) 



We can see this validated in Tables 1 and 2. All en- 
counters with a positive energy change show an increase in 
eccentricity, while encounters that reduce the energy of the 
binary's orbit also reduce its eccentricity, the magnitude of 
the change in eccentricity being approximately proportional 
to the disc mass for a given disc density profile. 
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Figure 18. A Jorbit / ^^orbit of the disc material at radius r plot- 
ted against r for the prograde (a — top) and retrograde (b — bot- 
tom) coplanar encounters. Points with A^'orbit > ^re shown as 
filled squares whilst those with A^'orbit < ^re shown as open 
squares. 



6 CONCLUSIONS 

We have used a reduced three-body method to evaluate 
the energy and angular momentum transfer between pro- 
tostellar disc material and the orbit of a perturbing mass. 
Such interactions are expected to be relevant to the early 
stages of binary formation and in the dynamical evolution 
of compact stellar systems. We have concentrated on the 
response of a disc around a central star to the parabolic 



fly-by of a perturbing mass (without disc). Evidence for the 
non-coplanarity of discs within binary systems (Hale 1994) 
prompted us to investigate several non-coplanar encounter 
orientations as well as the coplanar retrograde and prograde 
encounters. 

In almost all encounters we have found that both the 
energy and angular momentum transfers are dominated by 
disc material that becomes unbound from the system during 
the encounter. The only exception is in the case of the angu- 
lar momentum transferred by material at large radius in the 
prograde encounter, here dominated by material remaining 
bound to the primary - the contributions from the unbound 
and captured material approximately cancelled each other 
out. We have found no correlation between the original lo- 
cations of the material that yields peaks in energy transfer 
and that producing angular momentum peaks (Figs 4 and 
11). This suggests that a different mechanism is responsible 
for the two effects. Temporal analysis of the peaks in the en- 
ergy transfer showed that close encounters with the primary 
or secondary masses (or, occasionally, both) are responsible 
for peaks in energy transfer. Approximately one third of the 
energy transfer is in these peaks. Large transfers in angular 
momentum are caused by more distant encounters. 

In all cases where there is some prograde element to 
the motion, we find a corotation resonance region where 
the disc material loses energy and angular momentum to 
the binary orbit (Fig. 9). The maximum extent of this re- 
gion is out as far as r/fperi = 0.6 and thus for encounters 
more distant than this we expect the binary to become less 
strongly bound, and, if the original orbit were weak enough, 
the encounter could unbind the binary completely to pro- 
duce two single stars. The resonance has two contributors. 
The inner part is caused by the disturbance of disc material 
that remains bound to the primary whilst the outer part 
is due almost entirely to material that is captured by the 
perturber. This region clearly cannot be treated in a lin- 
ear fashion and the strength of the disturbance in the inner 
part implies that a linear treatment here may be doubtful. 
On the basis of this, we would expect linear calculations 
only to be valid out to a maximum radius of 0.25rperi (i.e. 
valid only for encounters with fperi > 4rdi3c) in the regime 
where M2 ~ Mi . We therefore conclude that the assump- 
tion of Ostriker (1994), that linear theory could be extended 
to fperi ~ 2rdi3c for equal mass encounters, is overly opti- 
mistic. When M2 <^ Ml we would expect this linear regime 
to extend to larger radii owing to the weaker nature of the 
interaction, but cannot quantify this effect with the present 
(equal mass) calculations. 

While CP93 reported that the prograde passage was 
most spatially destructive to the disc, this is not true ener- 
getically. We find that the retrograde encounter produces a 
larger energy and angular momentum transfer than the pro- 
grade case, by up to a factor of three or so. Also, the energy 
transferred during the encounter is greater than the binding 
energy of disc material exterior to periastron. The differ- 
ence is dependent on the orientation of the encounter, from 
a factor of two in the prograde case to a factor of five for 
the retrograde encounter. This implies that previous work 
(Clarke & Pringle 1991a, for example), in which the energy 
transferred is set equal to the binding energy of the disc ma- 
terial exterior to periastron, is somewhat pessimistic, so that 
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disc masses used in these calculations should be rescaled by 
a factor of between two and five. 

As a first step in analysing this problem, we have ig- 
nored the gas dynamics of the disc material on the basis 
of its limited effect during the fly-by. It is clear from our 
results, however, that not all of the stripped disc material 
ends up on diverging orbits. Outwards of r ~ 2rperi, about 
half of the energy transferred is by material undergoing close 
encounters with the stars. Such material is back-scattered, 
forming converging flows (Figs 6 and 13) in which we would 
expect shock-heating of the gas to occur with subsequent 
energy loss. Obviously this cannot be quantified with the 
current calculation, but requires careful analysis in future 
hydrodynamic calculations. We note that such close, con- 
vergent encounters are considerably less important inward 
of r ~ 2rperi and would therefore expect the current calcula- 
tions to provide the best estimates in the limits of encounters 
that are not highly penetrative. 

As is typical in such simulations, there is a multi- 
dimensional parameter space to be investigated, and such an 
exploration has to be an on-going project due to the massive 
amount of computing time required. It is our intention to 
begin this task, however, not only to establish some empir- 
ical relationships between the parameters that describe our 
models and the quantitative results, but also to determine 
the limits to which this reduced three-body approach may 
legitimately be applied before a more rigorous treatment is 
necessary. 

Up to now, we have ignored the reaction of the binary 
orbit to the disc material during the encounter, only cal- 
culating a bulk change at the end of the interaction. It is 
necessary to perform some simulations taking this into ac- 
count in order to evaluate how large the difference would be. 
Comparison of our post-encounter orbits with the parabolic 
initial orbit in the region near periastron where most of the 
interaction takes place suggests that any differences would 
be small. Viscous evolution of the disc material has also 
been ignored so far and an attempt to include this will also 
be made. 

It has been emphasized previously that systems at this 
stage of their evolution are deeply embedded in clouds of 
dust and gas which severely limit our ability to see them. 
Violent interactions such as those demonstrated here are an 
efficient method of removing/disturbing a large fraction of 
this material, thus dramatically affecting both the photom- 
etry of the object and its subsequent evolution. 
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